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Limited lithium supply is hindering the global transformation towards 


electrification and decarbonization. Current lithium mining can be 

energy, chemical and land intensive. Here we present an efficient and 
self-concentrating crystallization method for the selective extraction 

of lithium from both brine and seawater. The sequential and separable 
crystallization of cation species with different concentrations and 
solubilities was enabled by a twisted and slender 3D porous natural cellulose 
fibre structure via capillary and evaporative flows. The process exhibited 

an evaporation rate as high as 9.8 kg m°? h7, and it selectively concentrated 
lithium by orders of magnitude. The composition and spatial distribution 

of crystals were characterized, and a transport model deciphered the ion 
re-distribution process in situ. We also demonstrated system scalability viaa 
100-crystallizer array. 


Critical minerals, suchas lithium, have become a global development 
priority due to the exponential demand for such materials for batter- 
ies used in electric vehicles, consumer electronics and energy storage 
systems’ *. However, the sustainable supply of lithium can be a major 
barrier to support this energy transformation*”. Indeed, the price of 
lithium has skyrocketed by tenfold in the past decade, from US$6,000 
per ton in 2016 to nearly US$60,000 per ton in 2022, and the trend 
and fluctuation are expected to continue’. With the fast expansion of 
this market, the gap between demand and supply will keep increas- 
ing, asthe demand for lithium is expected to rise from approximately 
0.5 million metric tons of lithium carbonate equivalent (LCE) in 2021 
to more than 3 million metric tons in 2030 (refs. 7-9). It is therefore 
imperative to discover new approaches to boost new lithium produc- 
tion rate and capacity, not only from traditional salt lakes and continen- 
tal brines, which are accessible in very few countries, but also from new 
sources suchas low-concentration oilfield brines and seawater”. These 
new sources contain lower concentrations of lithium but can greatly 
expand lithium availability and ease geographical constraints’. For 
example, oceans contain ~5,000 times more Li than what is found on 


land, with -230 billion tons of Li available”, but the low concentration 
(-0.2 ppm) makes most extraction technologies unpractical". 
Currently, almost all lithium mining and production occurs in Latin 
America, Australia and China, and more than 50% of global lithium is 
currently produced from continental brines using large multi-stage 
evaporation ponds”. The world’s largest pond in Chile produced 20% 
of the global Li supply in 2019, but this evaporative process takes a 
long time (15-18 months) and requires a large land area (-32.3 km’) 
(ref. 16). Additionally, the operation is known to negatively impact 
local groundwater and ecosystems and has limited capability of pro- 
cessing brine sources with low lithium concentration (for example, 
<100 ppm) (refs. 17-20). Recently, new processes that perform direct 
lithium extraction have been developed by flowing brines through Li- 
selective materials using processes based on adsorption, ion exchange, 
electrodialysis and solvent extraction” ”™. These methods show poten- 
tial for increasing Li recovery and tapping into lower-concentration 
sources, but they are facing some major barriers in terms of selectivity, 
energy and chemical uses, and high cost. For example, the emerging 
electrodialysis with lithium-selective membranes relies ona selective 
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Fig. 1| Interfacial crystallization and spatial distribution of ion products on 
the twisted 3D fibre structure. a, Schematic of the 3D spatial crystallization 
process. At the bottom position (Bottom), mixed Li and Na ions are carried 
upward via capillary flow. At the middle position (Middle), NaCl saturates and 
crystallizes, while Li species continue to move upward with the water flow. At 
the top position (Top), Liions precipitate when the saline water completely 
evaporates. b,c, The Li/Na ratio increased along the 3D fibre crystallizer, with 

Li concentrated by >39 times at the top when using brine (b) and >675 times at 
the top when using seawater (c). Note that all Li/Na ratios along the height of the 
crystallizer were measured after removing the salt shell (details in the following 


section), exhibiting the Li separation performance within the crystallizer’s 
porous structure. d,e, A Na-rich salt shell forms on the outside of the crystallizer, 
which is easy to separate and harvest physically. d, Easy separation between 
cellulose fiber and the formed salt shell. e, Top-view of the separated salt shell. 

f, The Li/Na ratio is >200 times higher than that in the Na-rich salt shell, indicating 
the spatial separation of Li and Na in the radial direction. g, Comparison of the 
average water ER at the SQM Li production site in Salar de Atacama, Chile (22 °C, 
RH 22%), and the average ERs of a100-fibre-crystallizers array in the dark, under 
0.5 sun irradiation, or under 1.5 ms? wind. 


ion exchange mechanism, which can result in incomplete extraction 
of lithium. A study reported a recovery rate of -20% after 3 h, and the 
Li/Na selectivity was only 1.3 (refs. 26,27). In addition, electrodialysis 
employs an electrical field to drive ion exchange, which can be energy 
intensive. Solvent extraction is another method gaining attraction, as 
studies reported crown ethers have stronger complexes with Li* over 
other cations with a moderate Li/Na selectivity ranging from 2 to 5 
(ref. 28). However, current development is hindered by high material 
costs and the need for large amounts of energy and organic solvents”? =”. 
Moreover, after magnesium has been removed using chemical precipi- 
tation or electrocoagulation, Liseparation from monovalent ions (such 
as Na) becomes a major barrier” *°. 

In this Article, we report a passive method for fast and selective 
extraction of lithium from saline water based on twisted cellulose 
fibre crystallizers. By taking advantage of the high solubility and 
mobility of lithium salts in water”, we fabricated the crystallizers 
with unique porosity, tortuosity and hydrophilicity, which enabled 
selective extraction and spatial separation of lithium from other ions 


for simple harvesting (Fig. 1a). We show this crystallizer can accelerate 
the evaporation speed by more than 20 times compared with evapora- 
tion ponds, potentially reducing the operation time from >15 months 
to<20 days if properly scaled up. Moreover, compared with traditional 
evaporation, the process can selectively concentrate lithium by over 
39- and 675-fold from continental brine and seawater, respectively. 


Results 

The sequential and separable precipitation of lithium from sodium 
of different initial concentrations and solubilities is made possible by 
capillary and evaporative flows within the fibre crystallizer. The porous 
structure can raise water by -0.6 m via the capillary force and allows for 
fast water evaporation on the side surfaces. As water evaporates, salts 
with higher concentrations and lower solubilities, like NaCl, crystallize 
at lower heights of the fibre crystallizer, while salts with lower concen- 
trations and higher solubilities, like LiCl, move further upward and 
precipitate near the top (Fig. 1a). Therefore, with sufficient height of 
the twisted fibre structure, different salts reach their saturation points 
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at different locations and precipitate sequentially along the crystallizer, 
providing enough spatial separation for recovery. Figure 1b shows that, 
during evaporation operation, Li/Na wt. gradually increases along with 
the crystallizer height (H), and after 60 h, the Li/Na wt. (that is, 77%) at 
thetop of the crystallizer is 39 times higher than the typical continental 
Li-containing brine (that is, 2%). The average selectivity of Li versus Na 
onthe top 10% length of the crystallizer was 34.5 + 5.1, while the aver- 
age selectivity over the whole length was 7.1 + 1.2. In total, the lithium 
recovery by the whole length of the crystallizer was 87.2 + 1.9%, in which 
9.2 + 4.1% was accumulated on the top 10% length (Supplementary 
Fig. 1). A best recovery strategy that balances purity and quantity would 
consider the final product need, and multi-stage evaporation could be 
implemented to increase both. When using seawater, the crystallizer 
could raise Li/Na wt. by >675 times (that is, from 0.0015% to 1.02%), 
indicating that such spatial crystallization process can fit a variety of 
target ion concentration conditions (Fig. 1c). With a two-stage system, 
the final Li/Na wt. could reach a level comparable with brine-based 
lithium extraction (Supplementary Fig. 2). We provided more discus- 
sions in Supplementary Information regarding the theoretical and 
practical considerations (Supplementary Note 1). 

The spatial separation was also observed along the radial direction 
(R), with a Na-rich salt shell formed outside the crystallizer (Fig. 1d,e). 
Compared with a very high Na content in the salt shell, the Li/Na wt. 
was >200 times higher inside the crystallizer structure (Fig. 1f). This 
enables easy physical separation of Na salts from the crystallizer and 
further facilitates the recovery of concentrated Liin the porous struc- 
ture (Supplementary Movie 1). This spatial separation is a result ofion 
re-distribution within the porous material, which is influenced by the 
competition between advection and diffusion. With a high evapora- 
tion rate (ER) at the surfaces of the crystallizer, advection is dominant, 
and salts with higher concentration (NaCl) reach saturation first and 
crystallize at the air—solid interface where drying occurs. When the 
drying of the water expands from outside to the centre, amorphous 
Li salts start to precipitate in the centre. This spatial distribution of 
different salts allows the selective accumulation and separation of 
lithium from other cations. 

Inthe current design, the crystallizer features a very high height/ 
radius (H/R) ratio of ~600, which enables a large surface area and fast 
ER. A single device (projected area of 0.031 cm’) under an ambient 
laboratory condition showed an ER of 382 kg m° h”. An array of 100 
crystallizers (projected area 225 cm?) obtained an area normalized ER 
of 3.4 kg m°? h? (Supplementary Fig. 3). With 500 W m”? solar radiation, 
the ER improved to 3.7 kg m°? h”. When tested under the environmental 
condition simulating the actual weather in the Sociedad Quimica y 
Minera de Chile S.A. (SQM) LCE production site at Salar de Atacama, 
Chile (22 °C, relative humidity (RH) 22%, 1.5 m s” wind, dark), the 100 
crystallizers array showed an ER of 9.8 kg m7? h*, which is -27 times 
faster than the values reported by the largest evaporation ponds for 
lithium extraction (3,200 mm per year, SQM, Salar de Atacama, Chile)” 
(Fig. 1g). This suggests that the system can greatly shorten the opera- 
tion time of the evaporation ponds, potentially from >15 months to 
<20 days if properly scaled up. The increased ER may further improve 
ionselectivity per our mechanistic model, and details can be found in 
later sections. 

The demonstrated crystallizer is made from natural cellulose 
fibres fabricated with four twisted strings (Fig. 2a). A3D X-ray scanning 
image shows the cross-section of the 3D structure, where four strings 
are twisted tightly with an angle of ~23-24° to the vertical (Fig. 2b, Sup- 
plementary Fig. 4 and Supplementary Movie 2). When a 70-cm-long 
crystallizer was partially dipped in saline water, evaporation occurs 
continuously along the water lifting height (H) (Fig. 2c). Figure 2d 
shows the crystallization from brine after 60 hin ambient conditions. 
It was observed that most of the NaCl crystals formed surrounding the 
surface of the structure at around 1/3 height, while salt layers became 
gradually thinner upward (Fig. 2d). A time lapse video records that 


different original salt concentrations may result in different crystal- 
lization heights (Supplementary Movie 3). 

To decipher the composition and the mechanisms of sequential 
crystallization of different salt species, we characterized the salt den- 
sity distribution on two representative sections, one is at the very top 
(location 1) and another is at the top of the thickest salt accumulation 
(location 2) (Fig. 2d). The 3D X-ray images reveal the very different salt 
densities and distributions between the two locations (Fig. 2e). The 
most observable difference is in the longitudinal direction, the salt 
density was much higher at the location 2 than location 1, as evidenced 
by the colour change from red to blue. In addition, there is a density 
gradient in the radial direction, with lower densities inside the struc- 
ture, while higher densities outside the structure. This supports the 
hypothesis that Li salt is enriched inside the structure, while Na salt is 
enriched outside the surface. This distribution was further confirmed 
by other advanced elemental scanning methods. Overall, along the 
H direction, much higher Li/Na (>77%) was identified on the very top 
location, while lower Li/Na (-1.0-4.3%) was found in the lower 0.85 H 
section, with large NaCl cubic crystals dominating the formation ofa 
thick salt shell. 

We developed a mechanistic model to describe the ion transport 
inside the slender porous structure, which continuously draws liquid 
from the reservoir and evaporates along its length (Fig. 2g). The model 
predicts the concentration of Na (Cya) or Li (c,;) as a function of time t 
and the position zalong the water lifting height, H (Fig. 2h). The detailed 
model construction can be found in Methods. The predicted concen- 
trations critically depend on the effective Péclet number 
Pe = jH2/(pDR), a dimensionless parameter quantifying the ratio 
between advection and diffusion. Given the very high H/R, advection 
is dominant, causing an exponential increase in the concentration of 
both ions at the top of the crystallizer, c (z = H) = c%e?4/®). Such a 
sharp buildup of concentration, combined with the fact that Nais much 
closer to saturation in the bulk solution (cit > cle), leads to 
rapid precipitation of sodium. Consequently, after a characteristic 
time (pR/2/) In (cS /ch.,), the cy, remains constant at the top of the 
crystallizer, in equilibrium with the NaCl salt. This allows Li, which is 
comparatively further from saturation, to increase its relative concen- 
tration (Fig. 2i) as evaporation keeps driving the solution through the 
porous NaCl salt shell. Model results are consistent with our experi- 
mental observations (Fig. 1b,c and Supplementary Figs. 5 and 6), 
showing also a pronounced buildup of Li mostly at the top 10% of the 
wet crystallizer. Note that, when the ER increases, the maximum 
Li selectivity can greatly increase on the basis of the model prediction 
(Supplementary Fig. 7). This indicates that solar incidence and 
wind gust conditions in the field environment will accelerate 
evaporation and enhance efficacy compared with the lab condition 
(Supplementary Note 3). 

The experimental results support the hypothesis that in the radial 
direction (R), NaCl crystallizes outside the structure while LiCl con- 
centrates in the centre section (Fig. 3a). This is further supported by 
the cross-section scanning electron microscopy (SEM) with energy- 
dispersive X-ray analysis (EDX) scanning, where the Na elemental map 
shows a circular shape surrounding the crystallizer structure, andthe 
Li map shows a dispersed distribution across the crystallizer’s cross- 
section (Fig. 3b-d). When focusing on the area between the fibres and 
salt crystals (Fig. 3e, pink box in Fig. 3b), only salt crystals exhibit clear 
Na spectra, while the fibres showed no Na signal (Fig. 3f). In contrast, 
Li spectra were detected across the fibre areas, with the weakest signal 
at the bottom-right corner, where NaCl crystals are located (Fig. 3g). 
More images and discussion can be found in Supplementary Informa- 
tion, further confirming that NaCl crystals formed outside on the fibre, 
while LiCl formed amorphous salts across the fibre (Supplementary 
Figs. 8-10). Such a spatial distribution allows extraction and separation 
of Liand Na individually in practice, whichis demonstrated in crystal- 
lizer array testing (Fig. 5). 
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Fig. 2| Liis concentrated at the top of the fibre crystallizer. a,b, Zoomed 
optical image (a) and a 3D re-constructed image (b) show the twisted structure 
of the fibre crystallizer. c,d, The plain fibre crystallizer (c) and the salt-covered 
fibre crystallizer after Li extraction from the continental brine (d). Thick salt 
crystals are observed at the bottom 1/3 of the crystallizer, with a measured 
high concentration of Na, while thin salt layers distributed above that, with an 
increasing Li percentage towards the top. e,f, 3D density distribution at site1 
(e; pink-boxed site in d) and site 2 (f; blue-boxed site in d). Salt at site 1 (higher 


along the fibre) has a lower density, while salt at site 2 (lower along the fibre) has a 
much higher density, indicating that the lighter element, Li, distributed more at 
the top position while heavier element, Na, distributed more at the bottom. 

Note that the salt shell was removed before scan. g, Schematic of the geometry 
and parameters considered in the theoretical model. h, Prediction of the 

Liion concentration as a function of the position inside the crystallizer. 

i, Model prediction of the Li/Na ratio. 


We further characterized the salt crystal morphology to under- 
stand the dynamics of ion re-distribution and crystal formation. The 
SEM-EDxX investigations of a submillimeter-scale salt cluster can 
be found in Supplementary Figs. 11 and 12. Moreover, we employed 
electron energy loss spectroscopy (EELS) combined with S/TEM 
to investigate the crystal morphology and elemental distribution 


at the submicron scale. The S/TEM scanning images depict that 
square-shaped crystals are surrounded by amorphous precipitations 
(Fig. 3i). A clear boundary can be observed between the two types 
of solid in the zoom-in scan (Fig. 3j), indicating the presence of two 
different salt species. The EELS map confirmed this assumption of Na 
and Li distribution. On the other hand, Cl was present everywhere asa 
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Fig. 3 | Liis concentrated at the centre of the fibre crystallizer. a, Schematic 
showing that Li is concentrated inside the fibre crystallizer, surrounded by the 
Na-rich salt shell. b, SEM image of the cross-section. c, Na elemental map from 
the EDX, showing a circular shape surrounding the crystallizer structure. 

d, Li elemental map from the EDX, showing a dispensed distribution across the 
crystallizer’s cross-section. The Li elementary map was calculated by deducting 
the Na map from the Cl map (more details in Supplementary Note 2 and 
Supplementary Figs. 16 and 17)**"". e, SEM image of the area between the fibres 
and salt crystals, pink-boxed in b. f, Na elemental map from the EDX, showing 
that only the salt crystals emitted clear Na spectra, while the fibre area showed no 
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Na signal. g, Li elemental map from the EDX, showing that Li spectra were widely 
detected over the fibre areas, with the weakest signal at the bottom-right corner, 
where the NaCl crystals are located. h, Schematic showing that, in the centre 
Li-rich part, cubic NaCl crystals (blue) act as a core and LiCl form surrounding 
amorphous salt (pink). i, S/TEM image of the salt crystal sampled from the centre. 
j, TEM image of the blue-boxed area in i, showing a clear boundary between the 
two types of salt crystal. k, Layered EELS map confirmed two different elemental 
components aside the boundary. l-n, CI (1), Na (m) and Li (n) elemental map from 
EELS. Cl covered the whole scope, and Na distributed only at the right side of the 
boundary, while Li distributed mostly at the left side. 


counter ion (Fig. 3k-n). More images and discussions can be found in 
Supplementary Figs. 13-15. These results confirm that central square- 
shaped crystals are NaCl, while LiCl formed the amorphous precipita- 
tions in the peripheral region. The clear boundary between Na and Li 
suggests a sequential rather than a concurrent crystallization. Because 
NaCl has a lower solubility and higher concentration, NaCl crystallizes 
first, while LiCl is prone to precipitate around those NaCl cores gradu- 
ally when saline water dries out at the location. For the salt particles 
sampled from the salt shell, we did not observe any EELS spectra of Li. 
This confirms that most Li are localized inside the mesoporous fibre 
structure during sequential crystallization. 

The 3D twisted fibre crystallizer enables lithium separation by 
providing a high aspect ratio and a large surface area that allows staged 
crystallization of different salt species. The twisting structure fur- 
ther enhanced such features, as the spiral structure provides suitable 
porosities and pathways for saline water to move upwards, drying 
and crystallizing along the way. We fabricated the structure by firstly 
spinning the natural cellulose fibres into a tight string, then further 
twisted four strings together (Fig. 4a,b). For better comparison, cel- 
lulose fibres were also aligned to form a linear structure. Compared 
with the non-twisted control, the single-twisted string lifted water by 
102% higher, and the four-string twist crystallizer further increased 


the lifting by another 57% (Fig. 4c). Such twisting is believed to lead toa 
larger suction force to lift water higher, and we observed 351% increase 
in ER by the single twisted string compared with the untwisted string, 
and 108% increase in ER by the four-string twist (Fig. 4d). Our lab test 
also showed adding each pair of twisted strands may continue improv- 
ing water lift by ~10-15% (Supplementary Figs. 18 and 19). Detailed 
explanation can be found in Supplementary Note 4. Increasing the 
packing of fibres and strings by twisting would reduce their spacing 
and thus the effective pore size of the heterogeneous crystallizer. 
Since the permeability is directly proportional to the characteristic 
pore size of the medium, twisting could ultimately result in smaller 
permeabilities and larger values of the lifting height. Except in this 
section, all experiments were conducted using the four-string twisted 
structure as a single crystallizer. 

We also investigated the scalability of fibre crystallizers by 
constructing and testing a prototype array consisting of 100 fibre 
crystallizers (Fig. 5a). Such design can be easily modulated and scaled 
up for the quick deployment in various lithium production sites, 
which typically has the most favourable climate conditions for evapo- 
ration. Since NaCl would preferentially crystalize at the air—liquid 
interface, the formed Na-rich salt shells could be harvested by gentle 
vibration or scraping™. In contrast, LiCl stays within the crystallizer®. 
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Fig. 4 | Twist fibre structures achieve high water lifting and faster 
evaporation from brine water. a,b, Schematic (a) and photo (b) of the untwisted 
natural cellulose fibre (diameter 1.2 mm), one twisted string (diameter 1 mm) 

and four twisted strings (diameter 2 mm). c, Forward-looking infrared (FLIR 
image shows the water lifting heights of the untwisted string (18 + 1.6 cm), one 
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twisted string (37 + 1.5 cm) and four twisted strings (58 + 3.1 cm), indicating that 
twist structure can increase water lifting height. d, ERs of the untwisted string 
(24 +1.2 kg m° h”), one twisted string (108 + 15.2 kg m°? h”) and four twisted 
strings (225 + 29.7 kg m° h”), indicating that the twist structure can greatly 
facilitate water evaporation. 


Therefore, a wash-soak method (Fig. 5b) was implemented to fur- 
ther recover the Li-rich salt collected within the crystallizer. In this 
process, the outer Na-rich salt was washed off by quickly dipping the 
crystallizer (after outer salt shell was removed) in the water, and the 
inner Li-rich salt that retained on the crystallizer structure was then 
collected into a Li-rich concentrate by a thorough soaking. The Li/Na 
wt. in the recovered concentrate was found to be 55-105% higher than 
the pre-concentrated Li/Na wt. on the crystallizer structure, indicat- 
ing an improved overall selectivity of Li over Na (Supplementary 
Fig. 20). The estimated LiCl production from 100 fibre crystallizers 
is 2.7 kg per year after crystallization and wash-and-soak harvesting. 
Sucha process may save over 90% of the land use to reach a compa- 
rable yearly lithium production capacity of the SQM plant in Chile 
(-2.4 kt km” per year)), which has notable advantage for regions with 
limited space and environmental concerns”. A techno-economic 
analysis is needed to address the trade-off between the enlarged 
evaporation surfaces by stacking the crystallizers and the extra 
energy consumption for pumping brine higher, but it is beyond the 
scope of this study. 


Discussion 

The spatially separated crystallization enables effective Li extraction 
from saline waters without consuming large amounts of energy or 
chemicals, demonstrating a high efficiency, alow footprint and poten- 
tial process for scale-up. The process utilizes environmental energy 
as the driving force and can extract Li faster than the conventional 
evaporation pond method by an order of magnitude. The natural cel- 
lulose fibre-based crystallizer is low-cost and simple to fabricate and 
dispose, presenting an environmentally friendly method for lithium 
recovery. Further studies are needed to better understand the crystal- 
lization and separation mechanisms, and new materials and devices 
can be developed to improve separation selectivity and scalability. 
In this study we focused on the understanding of the ion crystalliza- 
tion process of LiCl via chemical analysis and in situ characterization, 
as for current evaporation pond operations, other ions such as Mg”* 
and Ca” are removed first before the final step of the most difficult 
monovalent ion separation. Further studies should investigate the 
potential interferences and removals of other ions present in natural 
waters to potentially simplify the current long trains of evaporation 
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Fig. 5 | Practical feasibility of the spatial crystallization system using 
simulated seawater. a, Photo of a prototype Li extraction array containing 

10 x 10 3D spatial crystallizers. Salt crystals are shown on blue cellulose fibre 
crystallizers. b, Lisalt recovery by the wash-soak method. After removing the salt 
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shell, the crystallizer was dipped into water to wash off and generate a Na-rich 
solution. Subsequently, the crystallizer was soaked in clean water bath for Li 
dissolution and recovery. 


process. We also hypothesize this approach could be used to extract 
other critical minerals such cobalt, nickel or uranium from brines and 
seawater. Overall, this new approach offers a sustainable pathway for 
securing the supply of the critical minerals and contributes to a circular 
low-carbon economy. 


Methods 

Interfacial crystallizer setup and ion analyses 

The base material used for fabricating the 3D spatial crystallizer was 
natural cellulose fibre yarn with a measured diameter of 2 mm. The 
cellulose fibre yarn was cut into 70 cm pieces and was then soaked in 
de-ionized (DI) water for 4 h to remove air bubbles in the pores to 
provide good liquid-solid contact along the whole structure. The fully 
wetted fibre crystallizer was then arranged vertically with the bottom 
~3 cm soaked in the saline water source with continuous replenishment 
during the experiment. After an extraction cycle, the crystallizer was 
cut into 2-4 cm pieces along the height and dried in a vacuum oven 
at 40 °C for over 24 h to fully remove the moisture. The easy-peeling 
salt shells were removed from the structure and dissolved in DI water 
for tests. Each individual crystallizer piece was also soaked in DI water 
separately for over 24 hto guarantee release ofall ions into the solution. 
The water samples were then diluted appropriately, and the ion con- 
centrations were tested by the inductively coupled plasma (ICP) after 
filtration by polyvinylidene fluoride syringe filters (25 mm, 0.45 um). 
Specifically, the Li* and Na* concentrations in all samples were meas- 
ured by the Thermo Scientific iCAP RQ ICP-MS. For system operated 
with synthetic brine, the samples were operated in the standard (STD) 
mode with 1,550 W plasma power and a carrier flow of 1.1 I min™ nebu- 
lizer gas. Samples with high salinity and low lithium concentration were 
measure in kinetic energy discrimination (KED) mode which uses a He 


collision Q-cell to remove interferences and allow low-mass ions such as 
Li to be detected accurately. To avoid errors caused by dilution, dupli- 
cate samples with different dilution ratios were tested for each sample. 


ER measurement 

The ER was calculated on the basis of the water loss of the reservoir. 
The reservoir containing DI water was placed on an electronic scale 
with real-time weight recording on a computer”. The whole system 
was placed in an environmental chamber (Controlled Environment 
Room, EGC) with controllable temperature and RH. The mass of the 
remaining water was recorded every 5 min, and the mass (y axis) was 
ploted against time (xaxis), where the slop of the plot was ER with a unit 
of g min”. The projected area of evaporators were also recorded. The 
ER for dark condition was tested under the ambient dark lab condition 
without sunlight exposure or wind assistance. The ER for 1 sun condi- 
tion was tested under the synthetic sunlight from a solar simulator 
with 300 W Xe-arc lamps (Newport), but without any wind assistance”. 
The ER for 1.5 m s” wind was tested under a measured wind speed of 
1.5 ms ‘that was generated from a commercial table fan, but without 
any sunlight exposure”. All the tests were conducted with an ambient 
temperature of -22 °C and an RH of 22%. For all the test, the blank refer- 
ence of evaporation from the reservior solely without any evaporator 
were also measured. The final ER for all the conditions were calculated 
by subtracting the blank reference. 


Advanced material characterizations 

The porosity of the crystallizer was measured by the water volume. 
The crystallizer was cut into sections with different lengths, where the 
weight was measured as mo. Then each section was immersed into DI 
water to fill all the pores. After removing the crystallizer from water, 
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the weight was measured as m,,... The absorbed water volume can be 
calculated by: 


Mwet — Mo 


Vi20 = P 


where p represents the density of H,O, which is 0.998 g cm” (20 °C). 
The volume of the crystallizer (V.,y<tatiizer) Was calculated by assuming 
the crystallizer as a cylinder. Then for different crystallizers with vary- 
ing lengths, we plotted Viystatiizer Ver SUS Vigo. The porosity is considered 
as the slope, which is 83.71%. 

Inthe micro- to millimetre scale, the morphology of the crystallizer 
before and after the Li extraction cycle was characterized by an environ- 
mental scanning electron microscope (SEM). Elementary maps for Na 
and Cl were acquired using the energy-dispersive X-ray spectroscopy 
(EDS) equipped in the same instrument”. However, the characteristic 
X-rays of Li are of too low energy (~55 eV) to be detected by EDS, so 
conventional elemental analysis/imaging methods can hardly detect 
Lisalts. Therefore, the elemental map of Li was adjusted by deducting 
the Na map from the Cl map, assuming that, if Na is absent, Li should 
be the counter ion pairing with Cl and balance the electric charge in 
salt crystals. To obtain accurate mineral separation results, we did not 
introduce other ions in these experiments (details in the next section). 

Inthe submicrometre scale, the characterization of the salt crys- 
tals was conducted by the EELS combining with the scanning/trans- 
mission electron microscopy (S/TEM). The sensitivity of EELS for Li 
is sufficiently high (0.2 wt.%) with the excitation of low energy loss 
appearing at 55 eV (Li K-edge). So the elementary map for Li, Na and 
Clinnanometres was directly acquired from the EELS, which revealed 
the clear distribution of Li and Na salts and supported the analysis of 
crystal formation. 


Salt species composition analysis 

Two synthetic saline water conditions were used for testing the Li 
extraction. Seawater condition was based on the typical ocean compo- 
sition: c(Li*) = 0.18 ppm and c(Na’*) = 10,800 ppm. Brine condition was 
based onthe typical geothermal brine composition at Salar de Atacama, 
Chile: c(Li*) = 1,500 ppm and c(Na*) = 76,000 ppm. 

After a Li extraction process from the brine condition, the salt 
density along the 3D spatial crystallizer was characterized by 3D X-ray 
images generated by Zeiss Xradia 520 Versa High-resolution 3D X-ray 
tomography Microscope. Two 2 cm sections were taken from both Li- 
rich and Na-rich region from the 3D spatial crystallizers for analysis, 
as boxed in Fig. 2d. Avizo software was used to reconstruct 3D volumes 
automatically or manually from 2D scan data, and used for dataintegra- 
tion, visualization and animation of the scan data. 


Lithium recovery and selectivity 
Lithium recovery was tested onthe basis ofa 10 x 10 crystallizer array, 
which hada projected area of 225 cm’. The gap between two crystallizer 
fibres was 1.5 cm. The test was conducted under the ambient tempera- 
ture of 22 °C, RH of 22% and wind speed of 1.5 m s7, and the solution 
employed was the typical brine composition at Salar de Atacama, Chile. 
During the test, the 10 x 10 array was placed in a container contain- 
ing 2,000 ml brine, where the initial and end mass and duration time 
were recorded. After the operation, three parallel crystallizers were 
collected from different array positions. Each collected crystallizer 
was immediately divided into four segments, from top to bottom, each 
with a length of 12 cm. Then the topmost segment was further divided 
into three subsegments, witha length of 3,3 and 6 cm. The salt shell was 
peeled off for each subsegment, and then the crystallizer and salt shell 
were fully dissolved in 25 ml DI H,O. The Li and Na concentrations were 
analysed for all the samples. The total Li mass was calculated over the 
whole crystallizer (with the salt shell). The Li recovery for each sector 
was calculated by dividing the Li mass in the crystallizer by the total 


Li mass (equation (1)). The Li/Na percentage was calculated for each 
crystallizer segment. At the same time, the selectivity was calculated by 
dividing the Li/Na percentage in the crystallizer by the Li/Na percentage 
in the initial brine solution (equation (2)). 


Li in crystallizer 


Total “100% W 


Li recovery = 


iig : 
—% in crystallizer 
a 


(2) 


Selectivity = ——W——. 
mm” in initial brine 


For different recovery lengths, the Li recovery was obtained by 
integrating the Li recovery over different segments of the whole crys- 
tallizer. While the Li/Na wt. for a certain length was calculated with the 
average Li/Na wt. among different sections. The mean and standard 
deviation for Li recovery and selectivity were calculated on the basis 
of the three parallel samples. 


Mechanistic model development 

A comprehensive model was developed for the transport of ions inside 
aslender, porous cylinder continuously drawing aliquid solution froma 
reservoir and evaporating along its length (Fig. 2g). Continuous evapo- 
ration along the length of a high aspect ratio structure provides high 
surface area for increased concentration enhancement from water 
loss. Within the slender crystallizer, the fluid flow is perpendicular 
to the evaporative flux, and the porous medium does not partially fill 
with air pockets due to the continuous liquid flux from the reservoir. As 
elaborated earlier, the theory accounts for the transport of ions through 
advection (which is induced by ambient evaporation) and diffusion, 
predicting the concentration of Na or Li ions asa function of time tand 
the position z along the water lifting height, H. The one-dimensional 
model considers the transport equation for a species of concentration c: 


att z OO =P an (4) 


where D is the diffusivity of the species through the porous cylinder 
and vis the vertical velocity induced by evaporation. An expression for 
v can be obtained from the Darcy equation for fluid flow in porous 
media, yielding 


U(Z) = Vo (1 = =) (5) 


The characteristic fluid velocity is vo = 2jH/ (pR), in which j is the 
ER, p is the solution density, and R is the radius of the crystallizer. 
Furthermore, the height H of the water column (Fig. 2g) is set by the 
competition between capillary forces, which promote the liquid rise, 
which is defined as H = y/@/ (ogv{k). The permeability k and porosity 
@ are intrinsic properties of the material of the crystallizer, while g and 
y are the gravitational acceleration and the surface tension between 
air and water, respectively. The 1D model assumes R/H < 1as in the 
experiments. 

The model equations are governed by a single dimensionless 
parameter, the Péclet number Pe = jH?/(oRD), which measures the 
relative predominance of advective and diffusive transport effects and 
has already been identified in other salt precipitation problems”. 
Additionally, the ratio of ion concentrations also depends on the 
saturation parameter cQ, /cit of sodium. The governing equations are 
solved numerically using an in-house finite-volume method witha 
two-stage implicit—explicit Runge-Kutta time integrator. Additionally, 
asymptotic analysis shows that, at high Péclet numbers Pe > 1, 
the leading order solution for the concentration of any ion is 
c(z,t)/c° = H/ (H — z) for O < z/H <1—e-2*/(R/) and c(z, £) (co = e% (PRI) 
for 1 — e~2/(PPRI) < z/H <1. 
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Data availability 


All data are presented in the article and its Supplementary Information. 


Source data are provided with this paper. 
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